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Doppler imaging of Speedy Mic using the VLT 

Fast spot evolution on a young K-dwarf star 
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Abstract. We study the short-term evolution of starspots on the ultrafast-rotating star HD197890 ("Speedy Mic" 
= BO Mic, K 0-2V, P ro t = 0.380 d) based on two Doppler images taken about 13 stellar rotations apart. Each 
image is based on spectra densely sampling a single stellar rotation. The images were reconstructed by our Doppler 
imaging code CLDI (Clean-like Doppler imaging) from line profiles extracted by spectrum deconvolution. Our 
Doppler images constructed from two independent wavelength ranges agree well on scales down to 10° on the stellar 
surface. In conjunction with nearly parallel V-band photometry our observations reveal a significant evolution 
of the spot pattern during as little as two stellar rotations. We suggest that such a fast spot evolution demands 
care when constructing Doppler images of highly active stars based on spectral time series extending over several 
stellar rotations. The fast intrinsic spot evolution on BO Mic impedes the determination of a surface differential 
rotation; in agreement with earlier results by other authors we determine an upper limit of \a\ < 0.004 ± 0.002. 
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1. Introduction 

Doppler imaging allows the reconstruction of surface maps 
of a sufficiently fast rotating star. To this end Doppler 
imaging makes use of the deformations passing through 
the spectral line profiles as the star rotates. The resolution 
of these maps depends on the spectral resolution, phase 
sampling and noise level of the spectra used as input. By 
making use of information modulated into the line profiles, 
Doppler imaging overcomes the diffraction limitations of 
direct and interferometric imaging techniques. 

Time series of Doppler images of the same object yield 
information on the evolution of surface features. This evo- 
lution is usually split up into "intrinsic" spot evolution 
(spot appearance and decay) and large-scale spot motion 
(tentatively related to meridional flows and differential ro- 
tation). Due to the considerable observational effort re- 
quired for constructing a high resolution Doppler image, 
in most cases only two images are available of the same 
object for direct comparison. As a result, the different spot 
evolution processes can usually not easily be disentangled. 
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As to intrinsic spot evolution, observational results 
on spot lifetimes are so far of a rather singula r nature , 
conce ntrating on a few well-observed objects ijHussainl . 
2002]); short-term studies, covering a few stellar rota- 
tions are rare. As an example, for the intensively stud- 
ied apparently single star AB Dor (K0V, P lo t = 0.51 d) 

only the pola r spot appe ars to survive time spans 

of years (e.g. iDonati et all ^999). We estimate other 
spots or spot groups, extending several dozen degrees 
on the surface, to persist longer than about 5 da ys 
l|Donati fc Collier Cameronlll997al:lDonati et allll999t as 
deduced from comparing their Doppler images of AB Dor) . 
Spots of approximately 10° size (close to the resolution of 
the available Doppler images) seemingly emerge and de- 
cay raisurnla^me^cales, i.e. over a few days. However, as 
e.g. lDonati et al.l l|l999|) emphasize, the phase coverage of 
the compared Doppler images must be checked carefully 
in order to derive significant lifetime estimates of individ- 
ual features. For the case of the rapidly rotating G dwarf 
star He 699 (fint. = 0.49 d) the available Doppler images 
l|Barnes et all Il998|) show that the large-scale distribu- 
tion of spots or spot groups remains stable over a month, 
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while features of approximately 10° size apparently evolve 
beyond recognition. 

Concerning large-scale surface flows, the presence of 
slow meridional flows in the solar convection zone (SCZ) is 
well established. According to helioseismic measurements, 
these flows extend at least several Mm deep into the SCZ 
l|Gizon et all 12003^ . However, long-lived sunspots do not 
show a regular meridional motion; their meridional dis- 
place ments are ap parently dominated by intrinsic evolu- 
tion llWohll 12002 V While there are possible indications 
of meridional mo vements of the spot patter n for a few 
evolved stars (e.g. lWeber fc StrassmeierL 12001^1 . at present 
no compelling observations of meridional flows exist for 
solar-like stars other than the Sun. 

As to differential rotation, the observational situa- 
tion is more encouraging. However, determining the ro- 
tation law of sunspots (i.e. their differential rotation) re- 
quires extended time series of observations. Due to the 
intrinsic evolution of sunspots it is impossible to derive 
their rotation law from two solar images alone. On the 
other hand, for suitable stars other than the Sun this 
appears to be possible due to the much larger number 
of spot features available in each surface image. The 
most prominent example is agai n AB Dor for which 
iDonati fc Collier Cameron! 1 1997aft measured its differen- 



tial rota tion strength, later c onfirmed by additional obser- 
vations ijDonati et all fl999t) . 

The angular velocity f2 of sunspot rotation as a func- 
tionof helipgraphic latitude B can be approximated by 
llBeckMl999ft 



sunspots(-B) « 2.9 • (l - 0.19 sin 2 B) ptx&d s _1 



(1) 



The precise coefficients of this rotation law depend on the 
size and lifetime of the sunspots considered, however, de- 
viations are at the level of a few percent. With the same 
accuracy, this law also describes the rotation of the solar 
near-surface plasma. 

Presently, a rotation law like Eq. ^ is also tenta- 
tively assumed for int erpreting observations of sta r s other 
than the Sun (e. g . IDonati k. Collier Camerorl Il997al 
iReiners fc Schmitil 12003ft . However, parameters to de- 
scribe the "strength" of differential rotation can be de- 
fined independently of the underlying rotation law. In the 
following we use 



(2) 



For the Sun, this definition results in a value of ag « 0.25. 
Whether a is a significant physical parameter to describe 
the dynamics of the outer convection zone is unclear at 
present. 

2. The target "Speedy Mic" 

Attention was drawn to BO Mic bv lBromage et ail ( 1992) , 
reporting a large flare observed during the EUV all-sky- 
survey of ROSAT; in retrospect it turned out to be the 



largest stellar flare observed during the whole mission. 
Their optical follow-up observations yielded a radial ve- 
locity of —6.5 ± 2.0 km/s with "no evidence for bina- 
rity" and an estimated v sini of 120 ± 20 km/s, earning 
BO Mic its nickname "Speedy Mic". A furth e r ana lysis 
of the same photometric data bv I Anders et"al] Jl993ft did 
not result in a well-defined rotat ion period; this issue was 
settled bv lCutispoto et ail 1 1997ft . who determined a value 
of 0.380 ± 0.004 days. 

In summary, its short rotation period makes BO Mic 
a promising target for Doppler imaging studies because a 
full stellar rotation can be observed during a single ob- 
serving night, allowing the reconstruction of the complete 
stellar surface. In this way the influence of intrinsic spot 
evolution during the observations required for a complete 
Doppler image can be minimized. 

A set of Doppler images, based on obse r vation s in July 
1998 has been published bv iBarnes et alJ l)200lft . As the 
authors state, the rather inhomogeneous SNR and phase 
coverage of their spectra did not allow detailed statements 
about the small-scale spot evolution or differential rota- 
tion on BO Mic. 



2.1. Stellar parameters 

2.1.1. Magnitudes and colours 

The unusually fast rotation of BO Mic for an apparently 
single dwarf star suggests a young evolutionary status. 
Actually, the photometric colours observed b y Cutispot o 
are not consistent w ith a main sequence sta r l|Cox[ EoOOft . 
The observations of ICutispoto et al ] ill 997ft yielded a vi- 
sual magnitude at maximum brightness of 



rfw « 9.32 ± 0.005 



(3) 



Together with the Hipparcos distance of 44.5±3.2 pc, im- 
plying a distance modulus of 3.24±0.16 mag, this results 
in an absolute visual magnitude of 

My w 6.1 ±0.1 . (4) 

The c olours of BO Mic determined by ICutispoto et alJ 
l)l997ft at brightness maximum are B — V = 0.92 ± 0.005 
and V - R = 0.56 ± 0.005. 

2.1.2. Evolutionary status 

Using the evolutionary models of ISiess et alJ |2000), in- 
cluding their colour-calibrations, the above values of My, 
B — V and V — R can be put into a consistent picture of the 
evolutionary status of BO Mic. Using models of 0.9±0.05 
solar masses we find 

Age 3.3±0.5 TO 7 yr 

Radius 0.9±0.05 R Q 

T oS 4750±50 K 

With T e ff denoting the effective temperature. The errors 
are determined from the evolutionary models based on 
the photometric errors given above. They should be con- 
sidered as rough estimates because of uncertainties of the 
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models and the fact that th e maxi mum brightness state 
observed bv lCutispoto et ah | lll997h does presumably not 
represent a completely spotless hemisphere of BO Mic. 

The above r adius estimate of 0.9 Rq for BO Mic re- 
sulting from the lSiess et alJ l)2000|) evolutionary models is 
roughly consistent with an estimate based on the observed 
rotation of BO Mic. An inclination angle of 70 ± 10° and a 
v sini of 134±10 km/s, both discussed in Sec. 03 result in 
an equatorial rotation velocity of w eg =143^^ km/s. This 
further corresponds to a radius of R =1.07 ^ Q14 Rq, using 
the rotation period of 0.380±0.004 days. In light of this, 
an inclination at the upper limit of the given inclination 
range (z = 60 — 80°) and a projected rotational velocity 
at the lower end of 124 — 144 km/s seem the be the most 
likely parameters. 

2.1.3. Lithium abbundance and kinematics 

Th e Li A6708 equivalent w idth of 220 ±50 mA determined 
bv iBromage et all l|l992h translates int o a Li abundance 
of 2.3± 0.4 dex, assuming T eff =4800 K. l|Soderblom et all 
Il993ch . Because of the large scatter of Li-equivalent widths 
and abundances for individual stars, this does not allow 
to closely constrain the age of BO Mic. As an example, 
the LiA6708 equivalent width ranges from about 100 to 
1000 mA at an effective tempe rature of 4800 K for th e 
Plejades cluster, age w 70Myr l|Soderblom et allll993c]) . 
For members of M34 (age w 200 Myr) an abundance of 
2.3 dex is approximately the upper limit detected at this 
effective temperature while for members of the Hyades 
cluster (age « 600 Myr) no lithium is significantly de- 
tected at effective temperatures below about 5400 K. In 
summary, the Li A6708 equivalent width of BO Mic indi- 
cates an age below "a few" 100 Myr which is compatible 
with the above value of 30 Myr, but does not help to con- 
str ain it any further. 

iMontes et al ] ll200lf) analyze the space velocity of 
a large sample of stars to check their membership of 
young stellar kinematic groups; their results concerning 
the membership of BO Mic in the Local Association 
(Pleiades moving group) are inconclusive, yielding no fur- 
ther hint of its age. 



3. Observations and data reduction 

3.1. Spectroscopy 

The spectral observations were performed at the VLT 
at the ESO Paranal using the spectrograph Uves 
(Ultraviolet and visual Echelle spectrograph). Uves was 
used in a dichroic mode, the spectral ranges thus covered 
were 3260 A to 4450 A and 4760 A to 6840 A for the blue 
and red arm, respectively. This permitted the observation 
of the ranges above 6000 A, well suited for Doppler imag- 
ing, and e.g. of the Can H&K lines, useful for diagnostics 
of chromospheric activity. We have not used wavelength 
ranges below 6000 A for the Doppler images presented 
here: The massively increased density of spectral lines for 



a K-star leads to severe line blending due to rotational 
broadening for an ultrafast rotator like BO Mic. 

The spectrograph slit width was chosen to correspond 
to 1" on the sky, offering a spectral resolution of A/ A A « 
40 000, adequate for the intended Doppler imaging. 

BO Mic was observed continuously from dusk to 
dawn during the two VLT nights (2002 August 2: 
JD 2452488.53 to JD 2452488.93, August 7: JD 2452493.51 
to JD 2452493.92) with the exception of a few reference 
star exposures. Only the spectra of Gl 472 were used as 
a template in conjunction with the presented Doppler im- 
ages. 

The exposure times of BO Mic were adjusted between 
140 s and 200 s depending on the observing conditions 
to avoid overexposing the CCDs. The raw exposures were 
summed in pairs to obtain the spectra finally used for 
the Doppler imaging. The resulting SNR for both nights 
ranged typically between 300 and 400 after rebinning to a 
binwidth of 0.08 A. 

The reduction of the spectra was carried out in the 
Iraf environment. After pre-processing the CCD images 
and removing scattered light, th e echel l e orde rs were ex- 
tracted by optimum extraction l|Hornel Il986fl . Optimum 
extraction is primarily intended for the extraction of 
"noisy" spectra, this would not be necessary in this case 
(even for the poorest spectra the SNR exceeded 200 after 
rebinning). However, its adaptive modelling of the order 
extraction apertures could deal efficiently with the spec- 
trum variations induced by the significantly varying seeing 
conditions (between 0.8"and 2.5"). 

One of Uves' CCDs (the MIT-CCD of the "red arm") 
shows strong variations (up to about 15%) in sensitivity 
on scales of a few pixels; these variations are visible on 
technical flatfields as a "brickwall pattern". In order to 
adequately correct for these small-scale sensitivity varia- 
tions it was crucial to individually extract a Hatfield for 
each spectrum using its own extraction aperture. 

3.2. Photometry 

The photometric observations were carried out by FvW 
using the 0.5m-telescope an d modular photometer at the 
Sutherland site of SAAO l|Kilkennv et all ll98Sl) . The 
standard SAAO reduction procedures were follo wed, de- 
scribed in the appendix to [Kilkenny et all l|l998h . 

BO Mic was observed continuously during 6 nights 
(2002 August 3-5, September 12, 13, 15), as far as 
weather conditions allowed. Unfortunately, the August 2 
and August 7 nights were lost due to bad weather at 
SAAO which inhibited almost simultaneous photome- 
try wit h our spectra l time series. Two E-rcgion stan- 
dards l)Menzies et all Il989j) . HD192844 and HD193132, 
were used as comparison star and check star, respec- 
tively. The estimated error of the differential photometry 
is 0.005 mag. 

The photometry did not reveal any flares during the 
observed time span. The September 12-15 lightcurve ex- 
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Fig. 1. Application of spectrum deconvolution (sLSD) to 
the "6400 A" wavelength region of two spectra of BO Mic. 
In the large panels the gray broad-lined curve represents 
the observed spectrum, the black spectrum shows the fit 
by sLSD. The narrow-lined spectrum is the template used 
for deconvolution (a synthetic Phoenix 5200 K spectrum 
without rotational broadening) . The inset panels show the 
resulting line profiles (broadening functions), their x-axes 
are annotated in units of 120 km/s. The lowermost panel 
shows the "continuum correction function" (CCF) used 
for both spectra; see text for details. 



hibits the same maxima and minima within ±0.01 mag as 
its August 3-5 counterpart fFig. lllfl . Due to the long time 
interval between the August and September photometry, 
corresponding to about 95 rotations of BO Mic, the latter 
did not add any new information related to the rotation 
period or the presented Doppler images. As a result, the 
September photometry of BO Mic is not further discussed 
in the following. 



4. Line profile extraction using sLSD 

Deconvolving spectra with the aim of extracting "aver- 
age line profi l es" fo r Doppler imaging was introduced by 
iDonati et al. | lll997h . Originally this was done in the con- 
text of magnetic Doppler imaging to allow the detection 
of the extremely weak signatures of po larization compo- 
nents. To this end IDonati et al ] lll997h apply a deconvo- 
lution to wide spectral ranges (e.g. spanning more than 
1000 A) with the aim of massively increasing the SNR of 




6120 6125 
Wovelength (A) 




6120 6125 
Wovelength (A) 



1.018 
1.000 — 
0.982 . 



Fig. 2. Application of sLSD to the "6120 A" wavelength 
region; the input spectra are the same as in Fig. Q ob- 
served on 2002 August 2 at rotation phases 0.337 and 
0.639, respectively. For this wavelength range an observed 
template spectrum (G1472) was used. 



the extracted line profiles compared to individual profiles 
from the input spectrum. 

Donati et al. coined the name "least-squares deconvo- 
lution (LSD)" for their implementation of this idea. We 
call our variant of the same idea "selective least-squares 
deconvolution (sLSD)". In contrast to the applications 
of Donati et al., sLSD concentrates on the analysis of 
comparatively narrow wavelength ranges, typically several 
10 A wide, containing only a handful of "strong" lines. 
While this comes at the price of requiring a higher SNR 
of the input spectra, sLSD allows to take the local char- 
acteristics of the used spectral range into account. 

4.1. The method 

The spectral line profiles of a sufficiently fast rotating star 
can be approximated by convolving the narrow-lined spec- 
trum of a slowly rotating star (called template spectrum) 
with a suitable broadening function. The broadening func- 
tion describes profile shaping mechanisms not confined to 
small regions of the stellar surface: rotation, limb darken- 
ing, macroturbulence and spots. Using this terminology, 
sLSD tries to approximate a broadening function common 
to the spectral lines of the select ed wavelength region . 

The description of LSD by Donati et al.l <)l997j) fo- 
cusses on a matrix, i.e. a linear formulation of the de- 
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convolution problem. Presently sLSD does not make use 
of the linearity of the proble m, but uses a Ley enberg- 
Marquardt algorithm instead ijPress et all Il992() . As in 
the case of Donati et al.'s LSD, the broadening functions 
of sLSD are discretely sampled and the value of each sam- 
pling bin is treated as one parameter to be optimized. It 
may appear as a disadvantage of sLSD's not-explicitly- 
linear treatment that it does not allow a simple error 
propagation and does not yield a formal error estimate 
as a byprodu ct, in contrast to a linear treatment (e.g. 
iBarnesl 12004). However, we think that this disadvantage 
is not very significant, because the errors of the broadening 
functions, which result from the deconvolution, are clearly 
dominated b y non-statistic al errors. This is supported by 
the results of lBarnesI i|2004h whose "deconvolved profiles" 
(corresponding to sLSD's broadening functions) show a 
"rippling effect" , most pronounced in the continuum re- 
gions surrounding the broadening functions, with an am- 
plitude massively exceeding their estimated statistical er- 
rors. Because of this "rippling effect" we have actually 
excluded these continuum regions around the broadening 
functions from the deconvolution. 

A truly significant error estimate of the broadening 
functions resulting from sLSD (or LSD) presently appears 
as an open question. We think that such an error estimate 
would require a systematic study partly based o n syn - 
thetic input data, e.g. al ong the l ines o f [Barnes! l|2004h . 
However, in addition to IBarnesl l|2004h . mismatches of 
the template spectrum presumably need to be considered: 
They increase the ill-posedness of the deconvolution prob- 
lem of sLSD (or LSD), leading to instabilities of the re- 
sulting optimization process which could be the reason of 
the "rippling effect" mentioned above. 

A different parametrization of the broadening func- 
tions, e.g. by Chebichev polynomials, can be easily imple- 
mented in sLSD. However, in practice this does not reduce 
the number of parameters needed to describe the broaden- 
ing function: The proper approximation of the large cur- 
vature at the transition into the surrounding continuum 
requires unwieldy high order polynomials. 



4.2. Application and ingredients 

For BO Mic we selected two wavelength regions for the 
extraction of line profiles by sLSD. They were inten- 
tionally chosen as quite different in character (Figs. ^ 
and EJ): While the "6120 A"-region is dominated by a 
single line (CaiA6122), the "6400 A" -region comprises a 
handful of stronger lines (the strongest are Fe I A6400 and 
CaiA6439). 

Figs. H an d d illustrate the application of sLSD to the 
two wavelength ranges. In the large panels the gray broad- 
lined curve represents the observed spectrum, the black 
spectrum shows the fit achieved by sLSD. The vertical 
dotted lines delimit margin regions contributing to the 
spectrum fit with reduced weights; the narrow flat regions 




6400 
Wavelength (A) 



1.11 
1-00 — 
0.89 



Fig. 3. Application of sLSD to a moderately wide spectral 
range around the "6400A" wavelength region. The input 
and template spectra are the same as in the lower panel 
of Fig.^ the latter is not shown here. The reasons for the 
rather poor fit to the spectrum are discussed in the text. 



around 6404 A and 6434A in Fig. ^ have been excluded 
from the fit (see below in Sec. I4.2.T|) . 

The narrow-lined spectra show the templates used for 
the deconvolution. In case of Fig. ^ a synthetic spec- 
trum generated for an effective temper ature of 5200 K by 
the s pectrum synthesis code Phoenix l)Hauschildt et all 
1999) has been used. The deconvolutions shown in Fig. |21 
make use of an observed template spectrum (Gl 472, 
K1V); spectra synthesized by Phoenix tentatively adopt- 
ing a wide range of element abundances and effective tem- 
peratures yielded significantly poorer spectral fits in this 
wavelength region. 

The application to relatively narrow wavelength inter- 
vals leads to several issues to be treated by sLSD: The 
treatment of end-effects at the borders of the wavelength 
interval, tentatively correcting for small template deficien- 
cies and stabilizing the convergence of the broadening 
function optimization. These issues are discussed below 
in Sees. EHH1 and 

4.2.1. End-effects and template deficiencies 

Fig. [3] shows an example application of sLSD using the 
same template as in Fig. ^ but yielding a poorer fit 
to the input spectrum. This rather poor spectral fit is 
due to mismatches with the template spectrum (leading 
to slightly wrong equivalent widths of some lines, most 
pronounced in the region around 6360 A) and regions 
containing echelle order overlaps not optimally merged 
(the region near 6460 A). To avoid the resulting defi- 
ciencies of the fit shown in Fig. three measures have 
been taken: (i) Narrowing down the wavelength region, 
carefully selecting the boundaries of the fitted wavelength 
range and excluding narrow intervals from the spectral 
fit (near 6403 A and 6435 A), (ii) Defining margin inter- 
vals at the boundary of the fitted range where the relative 
weights used for the ^-optimization of the spectral fit 
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are reduced; these margin intervals are marked by vertical 
lines in Figs.2] an dEl (iii) Introducing an adaptive "con- 
tinuum correction function" (CCF), shown in the lower 
narrow panels of Figs, ED 

The exclusion of narrow wavelength intervals from the 
spectral fit, mentioned in item (i), is motivated by tem- 
plate spectrum mismatches, i.e. apparently missing weak 
lines in these intervals. These missing lines in the syn- 
thetic template spectrum may be due to wrong element 
abundances supplied to Phoenix or wrong atomic data. 
However, we were unable to find element abundances lead- 
ing to an improved spectral fit, so wrong atomic data seem 
more plausible. Actually, the exclusion of the intervals is 
not necessary in this case, because it has only a marginal 
influence on the line profiles resulting from the deconvo- 
lution. 

Measure (ii), namely the margin intervals, would not 
be necessary if the selected wavelength range were sur- 
rounded by practically line-free regions, i.e. intervals of 
"undisturbed" continuum. However, such undisturbed re- 
gions are largely absent in the visible spectra of very fast 
rotating cool stars (which are the primary candidates for 
Doppler imaging). By introducing the margin intervals, 
sLSD reduces the influence of strong lines completely or 
partly outside the selected wavelength range. 

The most important measure to reduce this influence 
is the application of the above named CCF. The observed 
spectrum is divided by the CCF prior to the optimization 
of the broadening function. The CCF is chosen as a slowly 
varying function of wavelength, i.e. varying on larger scales 
than the individual line profiles of the spectrum. In this 
way the CCF adapts the continua of the observed and 
template spectra. This makes a continuum normalization 
of the input spectra prior to sLSD unnecessary. Such a 
normalization is usually very difficult for strongly rota- 
tionally broadened spectra. It remains so when using the 
CCF, but the process can be monitored during the appli- 
cation of sLSD, instead of being potentially "hidden" in 
the original spectra reduction. Additionally the CCF com- 
pensates for the effects of strong lines at the boundary of 
the wavelength region. This compensation is clearly visible 
in Figs. ^ an d El showing the deformations of the CCF at 
the boundaries. Finally the CCF tentatively corrects for 
some deficiencies of the template spectrum (wrong equiv- 
alent widths or missing weak lines). An example for such a 
correction can be seen when comparing the region around 
6410 A in Figs. □ and 

The CCF is iteratively constructed by alternating op- 
timization steps for the broadening function and the CCF. 
The CCF is chosen as a Legendre polynomial of fixed de- 
gree; this degree needs to be chosen sufficiently high to 
provide "flexibility" for the intended corrections. On the 
other hand, the degree should be adequately low to avoid 
the CCF adapting to deformations of individual line pro- 
files which are to be fitted by the broadening function. As 
an example, the degrees are 7 in Fig.Q] 5 in Fig. [3 and 3 
in Fig. 01 respectively. 




6120 
Wavelength (A) 




Fig. 4. Selected ranges from synthetic spectra, not ro- 
tationally broadened. The Phoenix 5200 K spectrum 
(black) has been used as the template spectrum for apply- 
ing sLSD to the "6400 A" wavelength range. The Phoenix 
4600 K spectrum (red/gray) is considered as a proxy for 
the spectrum of the spots on BO Mic, see text for details. 
For both panels, the quasi-continua have been normalized 
to unity, the red/gray spectrum is shifted in flux by -0.5. 



For the extraction of a time series of line profiles the 
CCF should be kept fixed, or its constancy closely moni- 
tored. Otherwise variations of the CCF, e.g. due to insta- 
bilities of the optimization could introduce "fake" varia- 
tions of the extracted line profiles. 

4.2.2. Stabilizing the convergence 

During the early phase of the broadening function opti- 
mization small-scale oscillatory solutions (typically with 
large bin-to-bin amplitudes) tend to appear. This ten- 
dency is further enhanced by the alternating optimization 
of the CCF and the broadening function. These oscillatory 
solutions massively disturb or even prevent the conver- 
gence of sLSD. To avoid this, sLSD employs a Tikhonov- 
regularization of the optimization of the broadening func- 
tion; i.e. it applies a "penalty function" to suppress strong 
small-scale gradients of the solution. The smoothing effect 
of the Tikhonov-regularization on the solution profile is 
adjusted by the weight of the penalty function during the 
optimization. This weight is fixed prior to the sLSD itera- 
tion; for input spectra of low or moderate noise it can be 
chosen very small. In this way it suppresses the oscillatory 
solutions without noticeably smoothing the final solution 
profiles, i.e. without a loss of potential starspot signatures 
in the profiles. 

4.3. Selecting the template spectrum 

The observed stellar spectrum results from a disk inte- 
gration and is in effect a weighted sum of contributions 
from undisturbed photospheric regions and spots. As a 
consequence, (at least) two intrinsic spectra take part in 
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the disk integration, but only one template is used for the 
deconvolution. Since the disk integration is dominated by 
the spectrum of the undisturbed photosphere, due to its 
greater brightness, the template spectrum for the decon- 
volution needs to be selected similar to it. In this way the 
intrinsic spectrum (or even different spectra) of the spots 
is not included as input information to the spectrum de- 
convolution. 

If the intrinsic spectra of the unspotted and spotted 
photosphere were identical, extracting a broadening func- 
tion common to the whole stellar surface would be strictly 
correct. Naturally, because of the strong temperature con- 
trast and other differing atmospheric parameters, these 
spectra are not identical. However, Doppler imaging based 
on spectrum deconvolution apparently is quite successful 
as demonstrated by the works of Donati et al. and also 
supported by the results of this paper. As discussed in 
the following, such a success can only be expected if the 
spot-photosphere temperature contrasts are either small 
or large, in a sense to be defined below. 

Fig- [3 illustrates a spot-photosphere temperature con- 
trast which is small compared to e.g. the Sun. The spec- 
trum of lower effective temperature is used here as a ten- 
tative proxy for the spectrum of the spots on BO Mic, 
the temperature contrast has been chosen as 600 K in 
this example. This approximately represents the situation 
assumed for the construction of the Doppler images dis- 
cussed in Sec. 15.31 because in the considered wavelength 
range the ratio of the quasi-continuum fluxes of the two 
spectra is about 0.5. 

In this case the dominating lines in the shown wave- 
length ranges are very similar for the spots and the undis- 
turbed photosphere. As a result, the line profile defor- 
mation by cold spots is predominantly due to the spot- 
induced change of continuum flux and approximately cor- 
rectly recovered by the deconvolution with a single tem- 
plate spectrum. Obviously, this does not apply to spec- 
tral lines strongly varying with temperature; consequently 
they should be avoided when using sLSD for the line pro- 
file extraction. 

If on the other hand a large spot-photosphere temper- 
ature contrast is assumed (e.g. 180 0K, as deduced fo r the 
highly active K-dwarf LQ Hya bv lSaar et aTll200lJ) . the 
flux emitted by the spots is very small compared to the 
surrounding photosphere (less than 10% in the example 
of LQ Hya, estimating from Phoenix synthetic spectra 
of corresponding effective temperatures). In this case the 
intrinsic spectrum of the spots has very little influence on 
the resulting broadening function. Again, the deconvolu- 
tion based on a single template spectrum should produce 
approximately correct results. 

5. Doppler imaging using CLDI 

5.1. The method 

We call the Doppler imaging algorithm used in the fol- 
lowing "CLEANlike Doppler imaging" (CLDI), it is based 



on the algorithm dev eloped b v IKiirsterl l|l993f) . CLDI is 
described in detail in IWolteil (|2004^) . further descriptions 
and studies of its behaviour will appear in a separate pub- 
lication (Wolter, 2005, in preparation). CLDFs rationale 
is very s imilar to the algorithm described and tested in 
IKiirsterl 111993ft . 

The mathematical formulation of the Doppler imag- 
ing problem resu lts in a matrix equation, D = R ■ I (e.g. 
IVoert et alJll987h . For this formulation, the line profiles 
are discretely sampled and assembled end-to-end giving a 
one-dimensional vector D. Likewise, the surface has been 
divided into a discrete grid of zones whose values have also 
been packed into a one-dimensional vector / . The matrix 
R is the so-called response matrix. At the given resolution 
it describes the observable line profiles as a function of the 
stellar surface. 

Doppler imaging is an inverse proble m ma king this 
equation ill-posed and ill-conditioned fe.g. lLucvll 1994ft . In 
a way, CLDI approaches the equation in a mathematically 
"daring" way: Since the inverse R _1 does not exist, it at- 
tempts a solution of the above equation by iteratively ap- 
plying the transpose of the response matrix R T to the ob- 
served line profile deformations. More precisely, it applies 
a slightly modified R T to the difference of the observed 
line profile deformations and their reconst ructed counter - 
parts at each iteration step. As discussed in lWolterl l)2004ft . 
this iterative "backprojection" of the line profile deforma- 
tions onto the reconstruction surface can be justified by a 
geometric interpretation of R T . 

CLDI is intended as an alternative to regularized op- 
timization schemes for the inversion of Doppler imaging 
problems, like maximum-entropy Doppler imaging. In con- 
trast to methods based on regularization, CLDI does not 
explicitly treat Doppler imaging as a continuous optimiza- 
tion problem. Instead it attempts a discrete approxima- 
tion of the observed line profile deformations in the way 
described above. In this way, inspired by the appearance 
of sunspots, CLDI models surface reconstructions of high 
contrast and discrete contrast steps. 

5.2. sLSD and CLDI 

The line profiles (more precisely the broadening functions, 
cf. Sec. I4.1fl extracted from average spectra of BO Mic us- 
ing sLSD show minor deviations from symmetry, an ex- 
ample is shown in Fig. |3J 

Phase- independent line profile asymmetries massively 
disturb the convergence of CLDI. The reason is that even 
weak deviations of this kind mislead CLDFs "backpro- 
jection" of the input line profiles used for approximately 
locating (groups of) spots on the surface. Such deviations 
are intrinsically impossible to fit for Doppler imaging, be- 
cause they lie outside the line center without migrating 
through the profile with rotation phase. We therefore have 
to make sure that the time series of line profiles is as far 
as possible free of such phase-independent asymmetries. 
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Fig. 5. Upper panel: Average line profile of BO Mic 
(stepped graph) extracted by sLSD from the "6400 A" 
wavelength region; the radial velocity is annotated in units 
of 120 km/s. The smooth graph shows a fitted analyti- 
cal rotation profile, its parameters are given in Table ^ 
Lower panel: Correction function £ computed from Eq.|6] 
See text for details. 

Such asymmetries demand care when interpreting the 
shape of the line pr o file in terms of stellar parameters (e.g. 
iReiners fc Schmitil Eool . For the purpose of Doppler 
imaging these asymmetries are less crucial, if they are 
small compared to the spot-induced line profile deforma- 
tions. 

In order to ensure the convergence of the Doppler 
imaging, the line profiles extracted from the spectra of 
BO Mic have been transformed to symmetric rotation pro- 
files in the following way, illustrated in Fig. [5] The shown 
profile was extracted from an average of all spectra of 
BO Mic taken on August 2, these spectra sample a whole 
stellar rotation quite evenly. In this way spot-induced line 
profile deformations should largely cancel out, assuming 
that the spot pattern does not change during the obser- 
vations. 

A symmetric rotation profile is determined by fitting 
an analytic rotation function / rot to an average line pro- 
file of the time series. / m t, is computed from a "rotation 
profile" G llOravMl992h 



G(v rad ) 



2(l-e)A3 + ±7reA 
nW(l - e/3) 



(5) 



by convolving it with a Gaussian of constant width. Here 
the symbol A = 1 — (v ra d/(c • W)) 2 has been used, com- 
puted from 



W = • A 



using the center wavelength of the considered line Ao- 

A correction function £ is computed from the average 
line profile b ave and the corresponding analytic fit / rot : 



£(«rad) = 



fmt{Vrad) 
frave(frad) 



(6) 



A sample correction funtion is shown in the lower panel 
of Fig. in that case the amplitude of £ is about 0.03, 



i.e. the correction amounts to at most 3% of the continum 
surrounding the line profile. Each individual line profile 
of the time series is then corrected by applying the same 
correction function 



6corr(«rad) = 6(^rad) ■ £(«rad) 



(7) 



The thus corrected 6 CO rr are use d as the input to CLDI. 
Using the same correction £ for the whole time series is a 
natural requirement; it avoids introducing artificial phase- 
dependent deformations of the line profiles, which could 
cause artefacts in the resulting Doppler images. 

Care must be taken that symmetric deviations between 
the average input profile 6 avo and f rot are not "corrected 
away" during this process. Such symmetric deviations are 
induced e.g. by polar spots (and other rotationally sym- 
metric spot configurations). As illustrated by Fig. such 
signatures were not detected for the average profiles of 
BO Mic. 

5.3. The images 

Two independent time series of line profiles were extracted 
from the spectra of BO Mic fSec. I4.2|l . the corresponding 
images are denoted by "6120 A" and "6400 A" in the fol- 
lowing. The images are shown in Fig. H3 the achieved line 
profile fits are shown for one of them in Fig. 

As stated above, the SNR of the input spectra ranged 
typically between 300 and 400, for some spectra it reached 
up to 500. Only for 3 of the spectra was the SNR below 
250. We have carried out tests of CLDI using synthetic 
input profiles with very similar reconstruction parameters 
as those used for our reconstructions of BO Mic (Wolter, 
2005, in preparation). These tests suggest a "critical" SNR 
range for the Doppler imaging reconstructions, mainly de- 
pending on the resolution aimed at and on the adopted 
atmospheric parameters. When the SNR of the input pro- 
files falls below this critical range, a serious deterioration 
of image details results. Above this level, the reconstruc- 
tion quality only increases weakly with increasing SNR. 
Our tests indicate that this critical SNR range lies between 
about 200 and 250 for the parameters adopted here. 

5.3.1. Image reconstruction parameters 

The image reconstruction parameters are summarized in 
Tabled The limb-darkening and v sinz values have been 
determined by fitting analytical rotation profiles / ro t to 
average line profiles (broadening functions) of BO Mic 
(Sec. 15. 2|) . It should be noted, that the errors of the fit 
parameters are presumably larger than the margins given 
in the table. This is mostly due to the asymmetries of 
the average line profiles, introducing uncertainties into the 
fit. 1 



Analyzing the same profiles in Fourier space 
iReiners fc Schmittl l2002ah yields v smi ~ 131 ± 2 kms 
However, since this method is also affected by problems of 
profile asymmetry, the error may well be larger. 



U. Wolter et al.: Doppler imaging of Speedy Mic using the VLT 



9 



Table 1. Parameters adopted for the construction of the 
Doppler images of BO Mic, e represents the linear limb 
darkening coefficient. The spot continuum flux is given 
relative to an undisturbed photospheric continuum flux of 
one. Note that the given margins designate the deviations 
between parameters adopted for different reconstructions, 
they are not proper error estimates. 



The adopted value of 0.5 for the spot continuum flux 
means that a surface element completely covered with 
spots (i.e. with a spot filling factor of one) emits half of the 
flux in the continuum around the synthesized line profile, 
compared to an uncovered surface element (spot filling 
factor zero). Adopting a lower value, i.e. a larger spot- 
photosphere contrast, leads to a slightly more "jagged" 
appearance of the spot groups of the resulting Doppler 
images, it has no significant influence on the reconstructed 
spot pattern. 

5.3.2. Inclination 

Adopting an inclination for Doppler imaging is commonly 
done by attempting reconstructions at different inclina- 
tions and selecting the one yielding the best line profile 
fits. Fig. 121 shows the line profile fit quality expressed as 
Xfi na | as a function of the adopted inclination for different 
datasets of BO Mic. These curves do not show a strong 
variation of xlnah this behaviour is not typical of CLDI, as 
we have tested for surfaces reconstructed from synthetic 
data. 

The minima of Xf\na\ f° r the August 7 reconstructions 
coincide at 70°. The August 2 reconstructions show less 
pronounced minima at lower values which do no agree. 
As discussed above, the radius of BO Mic tentatively de- 
duced from evolutionary models clearly suggest an incli- 
nation ^ 70°. As a result this value has been chosen for 
the reconstructions shown here. Actually varying the in- 
clination by as much as ±20° has only little influence on 
the resulting Doppler images. 

5.4. Surface resolution 

The resolution achieved by a Doppler image depends on 
the phase sampling, the spectral resolution and the noise 
of the spectra used for its reconstruction. Due to the one- 
dimensional projection of the surface structures onto each 
individual line profile, the resolution also depends on the 
reco nstructed latitude. Alth ough analytic estimates exist 
fe.g. l.Tankov fc FoingLll992Tl they are not much in use, be- 
cause the quantitative influence of the noise on the image 
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Fig. 8. Quality of line profile fits Xf\ n3 \ achieved by CLDI 
for different datasets of BO Mic as a function of the 
adopted inclination of the reconstruction star. The wave- 
length ranges and observation dates corresponding to each 
graph are annotated. 



resolution depends e.g. on the adopted atmospheric pa- 
rameters and requires simulations based on synthetic in- 
put data. In addition, potential errors of the line profile 
modelling, imperfectly treated line blends etc. can severely 
deteriorate the reconstruction quality. An extensive series 
of suc h simulations can be found in iRice fc Strassmeierl 
ll200ll) . We have studied the behaviour of CLDI based on 
synthetic input data llWolterll2004l 2005, in preparation). 
It has al ready been exte nsively studied for CLDI's prede- 
cessor in iKiirsterl l)l993|) with essentially the same results 
as for maximum entropy Doppler imaging. 

For the purpose of this paper, the resolution estimate is 
based on a comparison of the two presented sets of images. 
They are reconstructed from two disconnected spectral re- 
gions of rather different character (concerning the number 
of lines and the degree of line blending they contain). As a 
result, they can be considered as quite independent input 
data describing the same surface patterns. 

The images of Fig. E] belonging to the same input 
spectra (i.e. to the same dates) show corresponding fea- 
tures on scales down to about 10°xl0° on the surface; 
the correspondence becomes poorer within about 15° from 
the pole and the equator. This somewhat subjective esti- 
mate is substan t iated by the surface cross-correlation (e.g. 
iTonrv & Davisl Il979l) shown in Fig. El The correlation 
peaks are located within ±10° in the given latitude range. 
The cross-correlation was computed along longitude for 
each latitude strip of the compared Doppler maps. 



6. Results 

6.1. Observed spot evolution 

Comparing the August 2 and August 7 images, several 
spot groups are recognizable (Fig. EJ) : 
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Fig. 6. Doppler images of BO Mic, reconstructed by CLDI from the wavelength region given above each panel. The 
lower and upper map of each panel show the same surface, but rotated by 180 degrees. Black and dark gray/orange 
areas represent surface regions completely covered or 50% covered with spots, respectively. Unspotted regions are 
rendered as light shades. The subobserver longitudes of the observed phases are marked by short lines. A rotation 
phase <t> = corresponds to a subobserver longitude of ip — 0°; rotation proceeds with decreasing subobserver longitude. 
About 13 stellar rotations have passed between the images of the upper and the lower panels. 
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Fig. 7. Line profiles of BO Mic, extracted by spectrum deconvolution in the "6400 A" -region (red/gray) and fitted by 
CLDI (black). The reconstruction parameters are listed in Tab. Subsequent profiles of the time series are shifted, 
the flux scale is valid for the topmost profile of each panel. The rotation phases are given right of each profile, with the 
phase zero-point at JD^ = o = 2448000.05 Both observation nights were longer than the rotation period of 0.380 days, 
the hatched regions mark the cores of the profiles in the phase overlap interval of the beginning and the end of the 
night. 



Aq: A large irregularly shaped spotted area extending 
from about 120° to 240° in longitude and from about 
30° up to above 60° in latitude 

Ai: A smaller spot (group) centered slightly below 30° 
latitude and at about 290° longitude 

B: An extended area sparsely covered with spots, cen- 
tered close to 0° longitude and 30° latitude; its shape 
is reminiscent of a "A" 

Co: A spot group extending from about 50° to 120° in 
longitude and from close to the equator to about 60° 
in latitude. In the August 2 images this group is far 
more pronounced, showing a separation in two distinct, 
apparently quite coherent groups. 

C s : The small coherent southern spot only visible in the 
August 7 images at about 110° longitude and -25° lat- 



itude. It may be mirrored to the northern hemisphere 
in the August 2 images, appearing there as part of 
group C . 

The following evolution of the spot pattern appears to 
have taken place during the about 13 ( q ^dayso ) r °t a ~ 
tions between the August 2 and August 7 images: Group 
Ao undergoes rather weak changes but is shifted in lon- 
gitude by roughly 30°. Equivalently, Ao could be said to 
stay in place and fade at lower longitudes while growing 
at its high longitude side. Group Ai appears to stay in 
place, getting connected to Ao. Group B remains largely 
unchanged, preserving its A-shape and growing a weak 
"annex" to the north west at about 320° longitude and 
45° latitude. Finally, most of group Co disappears with 
possibly one component reappearing in the spot C 5 . 
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Fig. 9. Quantifying the Doppler image resolution: Cross- 
correlation of different-wavelength Doppler images recon- 
structed from the same input spectra (2004 August 2, i.e. 
the upper maps of Fig.^J). The left panel shows the cross- 
correlation of corresponding latitude strips of the com- 
pared maps; darker shades represent larger values. The 
relative maximum on each colatitudc (= 90°— latitude) is 
marked by a plus, the values of these maxima are plotted 
in the right panel. See text for details. 
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Fig. 10. Lightcurves computed from Doppler images of 
the "6400 A" wavelength region. The symbols represent 
the August 2 (o) and August 7 (x) reconstructions, respec- 
tively. The lightcurves are computed from the corresponding 
surfaces shown in Fig.|S| The annotation dates are JD-2440000. 
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Fig. 11. Lightcurve of BO Mic (2002 August 3-5), folded 
with a period of 0.380 days. Note that this folding pe- 
riod does not merge any portion of the data observed dur- 
ing different nights into a unique lightcurve. Annotation 
dates are JD-2440000: + = August 3, o = August 4 and 
□ = August 5. 



The apparent shift of group Ao towards increasing lon- 
gitudes (i.e. earlier rotation phases, keeping in mind that 
rotation takes place with decreasing sub-observer longi- 
tude), combined with the pronounced weakening of group 
Co at lower longitudes (i.e. later rotation phases) leads to a 
"darkening" at earlier phases and a "brightening" at later 
ones. As a result, the lightcurve minimum is shifted to ear- 
lier phases when going from August 2 to August 7; this 
can be seen in the synthetic lightcurves computed from 
the Doppler images, shown in Figure 1X771 

The available photometric lightcurve (2004 
August 3-5) could not be observed strictly parallel 
to the spectroscopic observations used for the Doppler 
imaging; even though it was observed on three consecutive 
nights, it contains unobserved gaps each lasting about 
1.5 rotations due to the short rotation period of BO Mic. 
However, the available lightcurve indicates that spot 
reconfigurations comparable to those seen in the Doppler 
images have occurred on time scales as short as about 
two rotations (namely between two consecutive nights). 
Due to the merely disk-integrated information about 
the spot pattern contained in the lightcurve, these spot 
reconfigurations cannot be further specified. 

The inferred fast changes of some of BO Mic's spots 
during the August 2002 observations on the one hand, and 
the snapshot-like information of the Doppler images (as 
well as the lightcurves) on the other, make definite state- 
ments about the spot lifetimes difficult. However, some 
inferences can be made: 

(a) Several spot groups of irregular shape, typically ex- 
tending 10-20° on the surface, located from equato- 
rial up to about 60° latitudes, reappear quite well- 
preserved after 13 rotations. 

(b) One quite coherent spot group of about 30° x 30° sur- 
face extension has substantially weakened after 13 ro- 
tations. 

(c) The only large spot group, quite coherent but of irreg- 
ular shape, extends more than 100° in longitude and 
about 50° in latitude. It reappears after 13 rotations 
with a comparable surface extension but significantly 
shifted and transformed on scales of up to 30°. 

A detailed comp arison between our D oppler images of BO 
Mic and those of lBarnes et alJ l|200lf) is difficult. The relia- 
bility of small-scale features in their images is presumably 
low, because of partly low SNR and partly rather inhomo- 
geneous phase coverage of their input spectra. However, 
the large-scale distribution of spots does show similari- 
ties. Barnes et al.'s 1998 images of BO Mic, as well as our 
2002 images presented here, both show spots over a wide 
range of latitudes including regions close to the equator. 
Both sets of images show one large spot group extending 
up to or very close to the poles, but neither the 1998 im- 
ages nor the 2002 images exhibit a polar spot. Our 2002 
images clearly show that the hemisphere opposite to the 
mentioned large spot group is considerably less densely 
spotted; this is supported by the modulation of the ob- 
served lightcurves rFier . fTTfl ■ Two differently spotted hemi- 
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Fig. 12. Cross-correlation of the August 2 and August 7 
Doppler images of BO Mic, reconstructed from the 
"6400 A" wavelength region (right column maps of Fig. [^J. 
The compared images have been smoothed by a two- 
dimensional boxcar extending 10°xl0° to account for 
their surface resolution estimated above. 



spheres also exist in some of Barnes et al.'s 1998 images. 
However, we think that due to the partly unfavourable 
phase sampling of the 1998 observations of BO Mic, no 
definite statement can be made about this. 

It must be kept in mind that all identified spot groups 
found in the Doppler images may comprise groups truly 
located on both stellar hemispheres, incorrectly recon- 
structed only on the northern hemisphere by the Doppler 
imaging. 

6.2. Differential rotation? 
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Fig. 13. Theory and measurements of differential ro- 
tation of (apparentl y) single dwarf stars. The c urves 
are adapted from iKitchatinov RiinVerl |l999) and 
extrapolated below a period of one day. Squares 
mark m easurements based o n comparisons of Doppler 
images dRice fc Strassmeieil [l996; B arnes et all [2000; 
iDonati et al.Lll999ll2003tl or the related meth od of "direct 
spot tracking" l)Colher Cameron et all l2002|) . The upper 
limit determined for BO Mic (Speedy Mic) in this work is 
rendered by the arrow's top. The arrowhead drawn outside 
the plot illustrates the possibility of anti-solar differential 
rotation (i.e. the equator rotating slower than the poles in 
terms of fi) which is not excluded by our observations. 



of 



a < 



-20 ± 10 deg 



0.004 ±0.002 



13.2 -360 deg 
i.e. at most 50 times smaller than the solar value. 



(8) 



The intermediate to high latitude spot group Ao is appar- 
ently shifted to larger longitudes (i.e. against the direc- 
tion of rotation) by about 30°, while features closer to the 
equator (including group Ai) roughly remain stationary; 
this result of a visual inspection is confirmed by the sur- 
face cross-correlation shown in Fig. ^] (see Sec. l5.4l for an 
explanation of the cross-correlation procedure). However, 
the surface shear visible in the cross-correlation map is 
not at all smoothly increasing from the equator to the 
poles as would be the case for a differential rotation law 
comparable to the solar case (Eq. . Instead, following 
the discussion of the preceeding section, we suggest that 
the apparent surface shear is predominantly due to the 
intrinsic evolution of the individual spot groups, i.e. their 
appearance and decay. 

As a result our 2004 Doppler images of BO Mic only 
allow us to estimate an upper limit of the differential ro- 
tation strength which could be veiled by the intrinsic spot 
evolution: Inspecting Fig. ^| yields a pole-equator shear 
estimate of Aip w — 20 ± 10° . The estimated scatter of 
±10° corresponds well with the resolution estimate of 
our Doppler images made in Sec. 15.41 This translates into 
a "relative strength" (cf . Eq. |2J of the differential rotation 



7. Summary and Conclusions 

We have acquired a time series of high-resolution spectra 
of BO Mic which is of outstanding quality regarding its 
low noise level and its dense and even sampling of rotation 
phases. BO Mic shows an exceptionally short rotation pe- 
riod which makes it a very interesting object for studying 
the spot evolution and possibly the differential rotation of 
ultrafast rotators (UFR). 

Our spectra allow the reconstruction of two Doppler 
images of BO Mic, seperated by about 13 stellar rota- 
tions. Each of these images is based on spectra observed 
during a single stellar rotation. This minimizes the influ- 
ence of intrinsic spot evolution on the images, which must 
be considered a possibility even during a single rotation. 

We have applied our spectrum deconvolution algo- 
rithm sLSD to extract line profiles from two separate 
wavelength ranges of the spectra. The Doppler images 
reconstructed from the resulting two independent input 
datasets agree on scales down to 10° x 10° on the stellar 
surface, excluding latitudes near the pole and the equa- 
tor. In this way we estimate the surface resolution of our 
images to about 10°. 
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While the large-scale spot pattern has survived the 
13 rotations covered by our images, several reconfigura- 
tions of spots have definitely taken place during that time 
span. These reconfigurations extend up to about 30°x30° 
on the surface. 

While our Doppler images show the timescale of these 
spot reconfigurations to be shorter than 13 rotations of 
BO Mic, our nearly simultaneous photometric observa- 
tions suggest that similar reconfigurations have even taken 
place between consecutive observing nights, i.e. during 
about 2.5 stellar rotations. 

Currently, for stars other than the Sun only very few 
observations of spot lifetimes on timescales of a few stellar 
rotations exist. By far the best-studied object concerning 
short-term spot evolution is the Sun. Compared to the 
solar case, "fast" changes of the spot pattern are not sur- 
prising: Pratically no sunspot group survives more than 
two solar rotations and definitely not without major re- 
configurations. 

The observed spectral line profiles of BO Mic contain 
remarkably rich information which agrees between the dif- 
ferent wavelength ranges. This fact and the resulting qual- 
ity of the Doppler images make our observations highly 
successful. However, one of our major aims could not be 
achieved, namely the observation of differential rotation 
on BO Mic and the determination of its strength. Given 
the large-scale stabilty of the spot pattern deduced from 
our Doppler images, the differential rotation can safely be 
said to be weak compared to the Sun. Unfortunately, fur- 
ther determination of the differential rotation strength is 
massively impaired by the above-discussed intrinsic evo- 
lution of the observed spots. 

Our determined upper limit of \a\ < 0.004 ± 0.002 for 
the strength of the differential rot ation agrees w e ll wit h 
a previously determined value by iBarnes et alJ l|200l|) . 
Combined with measurements of differential rotation on 
other fast rotating dwarf stars, shown in Fig. it con- 
firms the observed pronounced decline of differential rota- 
tion strength with decreasing rotation periods. As Fig. 1131 
also illust rates, this is well in accord wi th the rotational 
models of iKitchatinov fc Riidigerl (|l999lk so far support- 
ing the mean-field-based modelling approach for the con- 
vection zones of solar-like stars. 

Acknowledgements. U.W. acknowledges financial support from 
Deutsche Forschungsgemeinschaft, DFG - SCHM 1032/9-1. 
We thank D. Kilkenny for making the photometric observa- 
tions at SAAO possible and G. Cutispoto for kindly supplying 
photometric data of Speedy Mic. Finally, we thank our ref- 
eree J. Barnes for his careful reading of the manuscript and his 
concise suggestions which helped to improve it. 

References 

Anders G.J., Jeffries R.D., Kellet B.J., Coates D.W., 1993, 

MNRAS, 265, 941 
Barnes J.R., Collier Cameron, A., Unruh Y.C., Donati J.- 

F., Hussain G.A.J. , 1998, MNRAS 299, 904 



Barnes J.R., Collier Cameron, A., James D.J., Donati J.- 

F., 2000, MNRAS, 314, 162 
Barnes J.R., Collier Cameron, A., James D.J., Donati J.- 

F., 2001, MNRAS, 324, 231 
Barnes J.R., 2003, MNRAS, 348, 1295 
Beck J.G., 1999, Sol.Phys., 191, 47 

Bromage G.E., Jeffries R.D., Innis J.L., Matthews L., 
Anders G.J., Coates D.W., 1992, ASP Conf. Ser., Vol. 
26, 80 

Collier Cameron A., Donati J.F., Semel M., 2001, 

MNRAS, 330, 699 
Cox A. (ed.), 2000, Allen's Astrophysical Quantities, 

Springer, New York 
Cutispoto G., Kurster M., Pagano I., Rodono M., 1997, 

Information Bulletin on Variable Stars, 4419, 1 
Donati J.-F., Collier Cameron A., 1997, MNRAS, 291, 1 
Donati J.-F., Semel M., Carter B.D., Rees D.E. Collier 

Cameron A., 1997, MNRAS, 291, 658 
Donati J.-F., Collier Cameron A., Hussain G.A.J. , Semel 

M., 1999, MNRAS, 302, 437 
Donati J.-F., Collier Cameron A., Petit P., 2003, MNRAS, 

345, 1187 

Gizon L., Duvall Jr. T.L., Schou J., 2003, Nature, 421, 43 
Gray D.F., 1992, The observation and analysis of stellar 

photospheres, Cambridge Univ. Press, Cambridge 
Hauschildt P., Allard F., ApJ, 512, 377 
Home K., 1986, PASP 98, 1220 
Hussain G.A.J. , 2002, AN, 323, 349 
Jankov S., Foing B.H., 1992, A&A, 256, 533 
Kilkenny D., Balona L. A., Carter D. B., Ellis D. T., 

Woodhouse G. F. W., 1988, Mon. Notes Astron. Soc. 

S. Air., 47, 69 

Kilkenny D., van Wyk F., Roberts G., Marang F., Cooper 

D., 1998, MNRAS, 294, 93 
Kitchatinov L.L., Riidiger G., 1999, A&A, 344, 911 
Kurster M., 1993, A&A, 274, 851 

Lucy L.B., 1994, Reviews in modern astronomy, vol. 7, 31 
Menzies J. W., Cousins A. W. J., Banfield R. M., Laing 

J. D., 1989, SAAOC, 13, 1 
Montes D., Lopez-Santiago J., Galvez M. C, Fernandez- 

Figueroa M. J., De Castro E., Cornide, M., 2001, 

MNRAS, 328, 45 
Press W.H., Flannery B.P., Teukolsky S.A., Vetterling 

W.T., Numerical recipes in C (2nd ed.), Cambridge 

Univ. Press, Cambridge 
Reiners A., Schmitt J.H.M.M., 2002, A&A, 384, 155 
Reiners A., Schmitt J.H.M.M., 2003, A&A, 398, 647 
Rice J.B., Strassmeier K.G., 1996, A&A, 316, 164 
Rice J.B., Strassmeier K.G., 2001, A&A, 377, 264 
Saar S.H., Peterchef A., O'Neal D., Neff J.E., 2001, ASP 

Conf.Ser. 323, 1057 
Siess L., Dufour E., Forestini M., 2000, A&A, 358, 593 
Soderblom D.R., Jones B.F., Balachandran S., Stauffer 

J.R., Duncan D.K., Fedele S.B., Hudon J.D., 1993, 

ApJS, 85, 315 
Tonry J., Davis M., 1979, AJ, 84, 1511 
Vogt S.S., Penrod G.D., A.P. Hatzes, 1987, ApJ, 321, 496 
Vogt S.S., Hatzes A.P., Misch A. A., 1999, ApJSS, 121, 547 



U. Wolter et al.: Doppler imaging of Speedy Mic using the VLT 



15 



Weber M., Strassmeier K.G., 2001, A&A, 373, 974 

Wohl H., 2002, AN, 323, 329 

Wolter U., 2004, Ph.D. thesis, Hamburg 



